OPTIMAL DISTURBANCES IN BOUNDARY LAYERS SUBJECT TO 
STREAMWISE PRESSURE GRADIENT 


Anatoli Tumin 

The University of Arizona, Tucson, Arizona 85721 

and 

David E. Ashpis ' 

NASA Glenn Research Center at Lewis Field, Cle\’eland, Ohio 44135 


Introduction 

Laminar- turbulent transition in shear flows is still an enigma in the area of fluid 
mechanics. The conventional explanation of the phenomenon is based on the 
instability of the shear flow with respect to infinitesimal disturbances. The 
conventional hydrodynamic stability theory deals with the analysis of normal modes 
that might be unstable. The latter circumstance is accompanied by an exponential 
growth of the disturbances that might lead to laminar-turbulent transition. 
Nevertheless, in many cases, the transition scenario bypasses the exponential growth 
stage associated with the normal modes. This type of transition is called bypass 
transition. An understanding of the phenomenon has eluded us to this day. One 
possibility is that bypass transition is associated with so-called algebraic (non-modal) 
growth of disturbances in shear flows. ’ 
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A numerical analysis of spatial non- modal growth within the scope of the 
linearized boundary- layer equations for an incompressible flow over a flat plate was 
carried out in Refs. 3 and 4. Spatial analysis within the scope of the linearized Navier- 
Stokes equations ^uasi-parallel approximation of compressible and incompressible 
flows) was presented in Refs. 5-7. Recently, the method of Ref 4 was generalized for 

Q 

the case of compressible boundary layers. The main results of these theoretical 
models are as follows: 

• A system of counter-rotating streamwise vortices, which are periodic in the 
spanwise direction, provides the strongest growth of the disturbance. 

• There is an optimal spacing of the streamwise vortices, leading to the strongest 
effect. 

The effect of pressure gradients on the transient growth mechanism was considered 
within the scope of temporal theory by Corbett and Bottaro 9 and within the scope of 
spatial theory by Tumin and Reshotko. 7 Both studies were based on the quasi-parallel 
flow assumption. Tumin 10 analyzed the pressure- gradient effect for the Falkner-Skan 
profile within the scope of an analytical model when the spanwise wave number is 
very small. The pressure-gradient effect within the scope of spatial theory with 
nonparallel base flow and finite spanwise wave numbers has not been considered, yet. 

Another motivation for the present work stems from separation flow control on low- 
pressure turbines (LPTs). The performance of LPTs is strongly affected by the flow 
separation. There is a possibility of delaying the boundary- layer separation by tripping 
the boundary layer with the help of roughness elements or other devices. Usually, a trial- 
and-error method is used to determine an appropriate placement of the control elements. 
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This approach b time consuming and expensive. A recent investigation by Reshotko and 
Tumin 11 demonstrated that roughness- induced transition might be related to the transient 
growth mechanism. 

Periodically spaced in the spanwise direction, roughness elements generate a system 
of counter-rotating streamwise vortices. Due to a secondary instability mechanism, the 
streamwise vortices can lead to earlier transition to turbulence. They also provide a 
mixing enhancement due to redistribution of the streamwise momentum. Consequently, 
optimization of the streamwise vortices for maximum energy growth leads to 
maximization of the flow control effectiveness. In the present work, an analysis of the 
optimal disturbances/streamwise vortices associated with the transient growth mechanism 
is performed for boundary layers in the presence of a streamwise pressure gradient. The 
theory will provide the optimal spacing of the control elements in the spanwise direction 
and their placement in the streamwise direction. 

Governing Equations 

Because the flows of interest have relatively low Mach numbers, \\€ consider steady 
three-dimensional disturbances in an incompressible two-dimensional boundary layer. 
We choose the streamwise coordinate A' along the surface. The coordinate y will measure 

distance from the wall. We define a small parameter £ = ^jv / U re j-L re f that is the inverse 
square root of the Reynolds number, and v , ^ re j - , and L ye j- are viscosity, reference 

velocity, and reference length, respectively. The streamwise coordinate is scaled with 

Lr while the vertical coordinate v and spanwise coordinate z are scaled with 

i ej 
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^vZ re y / U re j- . The following scaling is assumed for the velocity disturbances 
u, v, and w , and the pressure p : 

u ~U v ~ eU r, w~eU p~e 2 pU f (1) 

rej ’ ref ’ ref ’ r r rej v J 

This scaling of the linearized Navier- Stokes equations and neglecting the curvature 
effects lead to the governing equations for Gortler instability, with the Gortler number 
equal to zero. We look for a periodic solution in the spanwise direction, with the 
corresponding wave number p . The governing equations for the amplitude functions can 
be written in dimensionless form as follows: 3 ' 4 



(2) 

( 3 ) 

( 4 ) 


-^-(Uw)+-^-(Vw)- P p = -~y~ P 


where U(x,y) and V(x,y) are the streamwise and normal velocity components of the 
base flow, respectively. The streamwise velocity, U(x,y), is scaled with [7^, and the 

normal velocity, V(x,y ) , is scaled with eU re j - . 

The following boundary conditio ns are applied to the solutions : 

y = 0 : a = v = w = 0 (6a) 
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y — » °o : 


u , w, p —> 0 


(6b) 


Equations (2)-(5) can be solved subject to boundary conditions (6a) and (6b) with 
prescribed initial velocity perturbations at x = Xq. 


Optimization of Energy Growth 

The authors of Refs. 3 and 4 employed an iterative procedure to find the optimal 
disturbances in terms of the maximum of the energy growth ratio G = E out / E jn , where 
E in and E out stand for the input and output energy norms. Andersson et aV used the 
same definitions of E in and E out as for the disturbance energy , 

Tmax . . 

E= J (u 1 +e 2 v 2 +£ 2 w 2 ^jdy (7) 

0 

whereas Luchini 4 employed the knowledge that the optimal disturbances are represented 
by streamwise vortices with corresponding output as streamwise velocity streaks, 

Anax 

E in= £2 \ (v 2 +W 1 )dy (8a) 

0 


liiaA 

E out= I u2 dy 


(8b) 


g = ^ou L = £ - 2 
E in 


Jmax 

1 

0 


lla A 

| u 2 dv 


3’max 


(8c) 


2 2 
v + w 


)dy 


As was shown in Ref. 3, the two definitions of the optimal disturbances lead to the 


same results at Reynolds numbers of 10 4 and higher. Because the iteration procedure 
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based on the optimization of ratio (8 c) provides significant simplification, we adopt it for 

2 

the following analysis. Because Eqs. (2)-(5) are independent of £ , the value of £ G is 
invariant with respect to the Reynolds number. 

Numerical Results 

Falkner-Skan Base Row 

We consider a Falkner-Skan family of boundary- layer profiles with free-stream 
velocity distribution U e = Cx m and corresponding Hartree parameter Ph = 2w/(w + l) . 
For the purpose of convenience, we have used the velocity scale U re j- =U eL = CL m and 
the length scale L re j- = L /(w + l) . The latter allowed the use of the conventional scaling 
of boundary layer solutions with H re y ~ -yj ^ T / ( w + 1) U 6 l — yJvL re j- / Uyef • 

Figure 1 shows the scaled energy ratio versus spanwise wave number ft for three 
Hartree parameters, P^ = -0.1, 0.0, and 0.1. The starting and the ending points, x jf1 / L 
and x out / L, are equal to 0.2 and 1.0, respectively. The Reynolds number Re^ in Fig. 1 
and what follows is defined as U eL L/v . One can see that an unfavorable pressure 
gradient ( p // < 0) leads to an increase in the energy growth while a favorable pressure 

gradient ( pfj > 0) leads to suppression of the transient growth mechanism. The latter is 
consistent with results obtained within the scope of parallel flow approximation 7 
Analysis of various starting points, x 7W / L , has shown that, in addition to an optimal 
spacing between perturbers, there is an optimal location from the leading edge (a similar 
result was observed in Ref. 8 for compressible boundary layers over a flat plate). 
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Example of LPT Conditions 


Volino 12 simulated low-pressure turbine (LPT) airfoil conditions in a low- speed wind 
tunnel. The test section was designed as a passage between two airfoils. The local free- 
stream velocity at a favorable pressure- gradient region was closely approximated by the 
following equation: 


Ue 

U exit 


1.48 


L c 


n0.214 


(9) 


where L s is the suction surface length and U exit is the nominal exit free- stream velocity 
based on the inviscid solution. The distribution (9) corresponds to a Falkner-Skan flow 
with the Hartree parameter pfj = 0.353 . 

Figure 2 demonstrates the energy ratio scaled with the Reynolds number Re^y = 
U ex j t L s /v versus the spanwise wave number scaled with Hi s =^JvL s / U exit . The 
ending point was prescribed at x out / L s = 0.444 while the starting points varied from 
0.11 1 to 0.289. One can see that there is an optimal starting point, x ir) / L The optimal 


velocity perturbation profiles at x out / L s = 0.444 (for //), x ir) / L s = 0. 1 1 1 (for v and w ), 
and pH i s =0.925 are shown in Fig. 3. 

The results indicate that we are dealing with a very strong favorable pressure gradient 
that suppresses the transient growth mechanism. For example, at a typical LPT cruise 
Reynolds number of 50,000, the transient growth will provide an energy amplification of 
less than 50. This is a relatively small number. If we take into account that, in practice, 
the perturber will not produce the optimal inflow field, the real amplification will be of an 
even smaller value. For example, in Blasius boundary layer, the theory predicts 
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amplification of 250 at the same Reynolds number of 50,000 (Ref. 9). Correlation 
between the transient growth factor and transition has not been established yet, therefore 
the effectiveness of the transient growth mechanism in preventing flow separation cannot 
be assessed quantitatively at the present time. 

There is a possibility of enhancing the transient growth mechanism by means of wall 
cooling. The effect of wall cooling was investigated by Tumin and Reshotko within the 
scope of a parallel flow approximation. In order to estimate possible increases of the 
energy ratio on a cold wall at a high favorable pressure gradient, we utilize the method of 
Ref. 7 for a compressible flow with local Mach number of 0.5 and Hartree parameter of 
0.353. The results are shown in Fig. 4. One can see that cooling of the wall might provide 
a tenfold increase in the energy ratio. 


Summary 

The Jesuits for the transient growth phenomenon within the scope of the linearized 
boundary- layer equations in the presence of a streamwise pressure gradient are consistent 
with previous results obtained within the scope of the parallel flow approximation and 
linearized Navier- Stokes equations. 7 A favorable pressure gradient decreases the noi> 
modal growth while an unfavorable pressure gradient leads to an increase of the 
amplification 

The example of a Falkner-Skan flow with a Hartree parameter = 0.353 
* 12 ’ 

corresponds to the experimental data “ and simulates the flow over a low-pressure turbine 
airfoil upstream of the separation point. At this pressure gradient, the transient growth 
mechanism is suppressed, and the energy amplification at low Reynolds number has a 
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small value. The theory of the transient growth mechanism predicts that it is possible to 
enhance the energy growth by means of wall cooling. The example within the scope of 
the parallel flow theory 7 demonstrates that cooling of the wall might provide a tenfold 
increase in the energy ratio. Future experiments on boundary-layer tripping accompanied 
by wall cooling will contribute to our understanding of the bypass transition mechanism. 

The method also predicts that there is an optimal spacing between perturbers and an 
optimal location from the leading edge. 

Consideration of the optimal velocity perturbations in Fig. 3 indicates that they are 
spreading across the boundary layer. This means that an array of generators localized on 
the wall will not provide excitation of the optimal disturbances. Therefore, the question 
of realizability of the optimal disturbances arises. For example, one can solve the 
receptivity problem for an array of generators on the wall and evaluate their shapes (or 
other parameters) to find the ones that provide disturbance profiles closest to the optimal 
ones. Another option is to solve the receptivity problem for distributed generators 
upstream of the starting point, x in , and to find the distribution of generators that leads to 

the optimal disturbances. The next option is to design a disturbance generator that 
directly affects the flow inside the boundary layer instead of perturbing the near-wall 
region only. These fundamental issues should be addressed in future research programs 
on the application of bypass transition mechanisms to separation flow control at low 
Reynolds numbers. 
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Fig. 2. Effects of the spanwise wave number p and the starting point xin I L on 
transient growth at the same conditions as the experiment in Ref. 12 ( p H= 0.353 , 
xout / Ls =0.444). 
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Fig. 3. The optimal velocity perturbation u at the ending point xout / Ls= 0.444 
and corresponding velocity profiles v and w at the starting point xin / Ls = 0. 1 1 1 . 
The parameters correspond to the experimental conditions in Ref. 12, (3 HLs = 0.925. 
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